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Abstract 



The Resistive Plate Chamber (RPC) muon subsystem contributes 
significantly to the formation of the trigger decision and reconstruc- 
tion of the muon trajectory parameters. Simulation of the RPC re- 
sponse is a crucial part of the entire CMS Monte Carlo software and 
directly influences the final physical results. An algorithm based on 
the parametrization of RPC efficiency, noise, cluster size and timing 
for every strip has been developed. Experimental data obtained from 
cosmic and proton-proton collisions at ^/s = 7 TeV have been used for 
determination of the parameters. A dedicated validation procedure 
has been developed. A good agreement between the simulated and 
experimental data has been achieved. 

Keywords: Resistive-plate chambers; Simulation methods and pro- 
grams 

1 The CMS Muon system 

The CMS (Compact Muon Solenoid) is one of the two general-purpose de- 
tectors located at the LHC (Large Hadron Collider). The main part of the 
detector is a superconducting solenoid providing axial magnetic field of 3.8 
T inside the coil and 1.8 T in the return yoke. The silicon tracker system, 
the crystal electromagnetic calorimeter, and the sampling brass/scintillator 
hadronic calorimeter are enclosed within the coil. The muon system is de- 
ployed outside of the solenoid and is situated between the layers of the steel 
return yoke. 

Three different technologies have been used to trigger and measure muons: 
drift tubes (DT) (pseudo-rapidity |r/| < 1.2); cathode strip chambers (CSC) 
(0.9 < \r]\ < 2.4) and resistive plate chambers (RPC) (|?7| < 1.6). The muon 
system consists of four stations in the barrel and in the endcaps. In the 
barrel part, each station consists of 12 DT layers packed on both sides by 
two RPC in the first and the second stations and one RPC in the third. 
The outermost DT station has eight layers and one RPC. In the endcaps, 
each muon station consists of six CSC detection planes complemented by 
the RPC in the three inner stations. A detailed description of the CMS 
detector can be found elsewhere [T]. 

The RPC system consists of 480 barrel chambers and 432 end-cap cham- 
bers and covers an active area of 2953 m^ f2]. This makes the RPC system 
the largest sub-detector at CMS. 

The ability to work at a high rate of ionizing particles up to 1 kHz/cm^ 
has been ensured by Bakelite double-gap RPCs operating in avalanche mode. 
A double-gap encompasses two single gaps mounted one on the top of each 
other with common readout strips situated in between them. The rate 
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limitations are due to the time needed to re-establish the electric field after 
an avalanche in a gas gap. The voltage drop is proportional to the average 
avalanche charge in the gas gap and the electrode resistivity. The analytical 
description of the rate effects in RPC is based on simple DC model (see [3], 
[Ij J [5] ) . The advantages of the double gap design are clear within the model, 
i.e. doubling the induced signal while keeping the charge released per gap 
at the same level than in a single gap. The electrodes are made out of a 
Bakelite which has lower resistivity than glass and the Bakelite resistivity is 
controlled by adding water vapours to the gas mixture [6]. 

Each chamber consists of two or three adjacent double-gaps (Figure [T]) 
held together by aluminium profiles. There are 1020 double-gaps in the 
barrel and 1296 in the end-caps. Each double-gap has up to 96 read-out 
strips. There are more than 10^ read-out channels in total for the whole 
RPC subsystem. The chambers in the barrel have rectangular geometry and 
the width of the strips varies from 2.3 cm for the innermost muon station 
to 4.1 cm for the outermost muon station. Chambers in the end-caps have 
trapezoidal geometry and the width of the strips varies from 1.7 cm at the 
lowest radii to roughly 3.6 cm at the highest radii. 

The RPCs have been used as dedicated muon trigger detectors. Due to 
the time resolution of about 1.8 ns, proper association between the RPC 
trigger information and the LHC bunch crossing has been ensured [6]. 
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Figure 1: Schematic layout of a RPC chamber composed of two double gaps. 



2 The CMS RPC simulation 

CMS experiment should be able to distinguish the detector related effects 
from the real physical processes. The CMS experiment studies different 
physical processes, detecting the particles originating from proton-proton 
or lead-lead collisions. A crucial part of the data analysis is the successful 
factorization of the detector related effects from the real physical processes, 
resulting in the collisions. For example one needs to know very well the 
geometrical acceptance and efficiency of the detector and the trigger system 
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in order to calculate the cross section for a given physical process. The CMS 
detector performance has been very well studied in details and has been im- 
plemented in the CMS simulation, thus simulation has been used to account 
for detector related effects, during data analysis. A dedicated Monte Carlo 
simulation of the CMS detector is a crucial part of the data analysis and 
the simulation result enters in the final physical result of the CMS collabo- 
ration. The correct simulation of the detector response is used also for the 
study of the upgrade of the CMS RPC system. In the past decades several 
RPC simulation techniques are proposed by different authors [8], [9], [10] . 

|12j . Most of them are focused on the simulation of the RPC properties 
,,ab initio", i.e. obtaining detector characteristics (charge spectra, efficiency, 
timing, etc.) from detector parameters (e.g. Townsend coefficient, gas gap 
width, number of gaps, etc.). These techniques are very important and use- 
ful for detector R&D. They can be used to simulate different RPC types and 
to optimize detector design in terms of gap width, gas composition, strip 
parameters. Even an approach for modeling of the RPC behavior, based 
on artificial neural networks is developed [13]. The method predicts the be- 
havior of the anode charge as a function of the environmental temperature 
and pressure. All of these techniques are good for describing the parameters 
of a given RPC type. Our task is to simulate 2316 individual detectors, 
all together with the detector electronics. All of these detectors are of one 
type in terms of Townsend coefficients, efficiency and front-end electronics. 
The detectors are working for almost a decade and the important detector 
parameters have been measured and well understood. To take advantage of 
our knowledge we adopted data driven parametrized simulation model. Col- 
lision data and cosmic data collected in 2011 have been used for estimation 
of the RPC efficiency and noise simulation parameters. The time resolution 
and cluster size distribution have been estimated from experimental data 
collected on cosmic rays [6]. 

The following parameters, taken from a dedicated database have been 
used in the default digitization algorithm: 

• efficiency for each readout strip; 

• cluster size distributions; 

• detector time resolution. 

• noise rate for each readout strip; 

In order to tune the MC simulation to the experimental data it is very 
important to update regularly the parameters used for the RPC response 
simulation. Presently, three sets of parameter values with different efficiency 
and noise rate parameters are used: 
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• default values : RPC efficiency 95%; strip noise 0.05 Hz/cm^; time res- 
olution 1.42 ns; cluster size distribution estimated from experimental 
data collected by cosmic ray telescope |7|; 

• parameters estimated using cosmic-ray data; 

• parameters estimated from proton-proton collisions data. 

2.1 The simulation chain 

CMSSW is a dedicated software framework for CMS event reconstruction, 
Monte Carlo (MC) simulation and data analysis. Simulation of the RPC 
response is an integral part of the entire CMSSW simulation chain which is 
used for the physical CMS analysis. The simulation chain in CMSSW has 
been organized as follows: 

• event generation - events have been generated by physics events gen- 
erators, interfaced by CMSSW. In the particular case Pythia event 
generator [14] has been used in order to generate Drell-Yan events 
with two muons in the final state [15]; 

• event simulation - CMSSW interface to GEANT4[16j has been used to 
simulate the passage of particles through the detector. The whole CMS 
detector has been described precisely in the simulation, taking into 
account the mechanical design, material budget and magnetic field. 
The RPCs have been described by their components i.e. aluminium 
housing, Bakelite plates, active gas volume, filled with the ,, standard" 
gas mixture [6] and readout strips. In the particular case of RPC 
simulation, the result of the simulation step is a simulated hit in the 
active detector volume; 

• digitization - in response to each simulated hit, one or few adjacent 
strips have been fired with a given probability. Timing information has 
been used to assign the detector response to a given bunch crossing. 
The intrinsic detector noise has been simulated on the same step; 

• trigger - Level 1(L1) trigger has been emulated and High Level Trig- 
ger(HLT) algorithms have been executed on this step; 

• event reconstruction - the reconstruction algorithms have been per- 
formed. 

2.2 The simulation parameters 

• Efficiency - the RPC efficiency has been estimated employing the seg- 
ment extrapolation method [2j. DT/CSC reconstructed track segments 
have been extrapolated to the RPC strip planes and used to look for 
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RPC hits in the surroundings of the extrapolated point. Efficiency has 
been computed as the ratio between the number of observed RPC hits 
and the expected number of hits from segment extrapolation. The av- 
erage chamber efficiency computed with this method has been assigned 
to all the strips of the same chamber. These values obtained from 
experimental data have been used for the RPC efficiency parametriza- 
tion. During the simulation, the decision whether the strip should be 
ffied or not has been taken according to the assigned strip efficiency. 

• Cluster size - the number of simultaneously fired adjacent strips (clus- 
ter size) has been calculated using an empirical cluster size distribu- 
tion for each chamber. The RPC intrinsic efficiency and cluster size 
are assumed to be independent from the muon momentum in our sim- 
ulation model. For our studies only high momenta muons (above few 
tens of GeV) have been taken into account. In a good approximation 
they behave as minimum ionization particles in our detectors |17j . For 
even higher momenta muons (in TeV range) the radiative effects are 
treated on the simulation step by GEANT. According to the experi- 
mental data [6], the clusters size has very weak dependence from the 
muon impact angle for angles smaller then 45° and thus that depen- 
dence is not implemented in the simulation. However, detailed analysis 
of the experimental data shows that the cluster size depends on the 
muon impact point on the strip. The dependence is parametrized in 
the simulations and the parameters are obtained from the experimen- 
tal data. The strip has been formally subdivided in five slices along 
its width, and the impact point of the muon has been used to find 
the fired slice coordinates. For each slice we have a different cluster 
size distribution and the random number, corresponding to the cluster 
size, is generated from the distribution. All the distributions have been 
obtained from the experimental data, using an algorithm explained in 
[18j . An example of such a distribution is presented on fig. |4j 

• Timing - the read-out of the CMS detector is performed every 25 ns. 
The true signals correspond to the signals within a 25 time window, 
synchronized with the LHC bunch crossing separation, so called central 
bunch crossing. The information from six consecutive time windows, 
two before and three after the central bunch crossing, has been col- 
lected to form an event. For proper bunch crossing assignment the 
time of flight of the particle, the signal propagation speed on the strip 
and the delay of the signal along the cables for each chamber have been 
taken into account. The time of flight depends on the momentum of 
the muon and has been calculated for each simulated particle. The 
signal propagation time along the strip depends on the hit position 
on the strip. It has been simulated assuming the signal propagation 
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speed equals to 66% of the speed of light. The value assumed for the 
signal propagation speed is consistent with the average shift of the sig- 
nal arrival time observed in data samples taken at different positions 
along the read-out strips for a set of chambers [19]. The uncertainty 
of the signal propagation speed is not an issue in the simulation, be- 
cause the time is measured in bunch-crossings, i.e. time quanta is 25 
ns. The intrinsic time resolution of the chambers altogether with the 
attached front-end electronics is less than 2 ns for all the chambers 
and the time resolution has been simulated as a Gaussian with a mean 
value ns and sigma 1.42 ns. The different channel-to-channel delays 
in a chamber are far below this value. A constant time offset for every 
chamber is added to account for different cable lengths, muon time of 
flight etc. The RPC signal synchronization is described in |20| . 

• Noise ~ there are two possible types of parametrization used for in- 
trinsic noise simulation. A parametrization where all the strips have 
a default value of 0.05 Hz/cm^ was used in the first half of the 2011. 
A realistic noise parametrization has been used afterwards. The in- 
trinsic RPC noise has been measured during cosmic runs |[6| and used 
to model the MC response. The estimation of the simulated noise 
rate for a given strip in a certain event has been based on the Poisson 
distribution with mean v: 

V = NjSstrt 

where Nj is the measured strip noise, Sgtr is the strip area and the total 
simulated time is t = ni,x.25ns, where is the number of simulated 
bunch crossings. 

The simulated noise signals are uniformly distributed in time windows 
around the central bunch crossing, including the central one. 

3 Validation 

The input to the digitization is efficiency, cluster size distribution, time reso- 
lution and noise rate, normalized by the area. The output of the digitization 
is a digi. The digi in the RPC case is a pair of two numbers. The first num- 
ber correspond to the read out channel of the front-end electronics connected 
to the fired strip and the second number corresponds to the time when the 
strip is fired, measured in bunch-crossing units (time quanta is 25 ns). The 
digi is used in the CMSSW framework for the consecutive Level 1 trigger 
emulation and on the reconstruction step. It is difficult to compare digis to 
the input parameters and in order to validate the whole simulation and the 
digitization in particular a dedicated analysis has been performed. In the 
analysis the output of the reconstruction of the simulated events has been 
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used to estimate efficiency from MC , cluster size distribution from MC and 
timing from MC. The results have been compared to the input parameters 
or to the real experimental date in order to estimate the goodness of the 
MC simulation and correctness of the MC model. 



3.1 Validation of the goodness of the simulation 

In order to validate the goodness of the simulation, Drell - Yan events with 
two muons in the final state with invariant mass greater then 60 GeV/c^ 
have been simulated using the Pythia event generator [14j. The generated 
events have passed the full simulation chain. A dedicated selection procedure 
choosing muons from Z decays has been applied both on the MC and collision 
data. For the purposes of validation, the efficiency of a given RPC has 
been estimated as the ratio between the number of events with at least one 
ffied strip corresponding to a simulated muon hit and the total number of 
simulated muons crossing the chamber. The comparison between the input 
efficiency parameters estimated from experimental data and the simulated 
efficiency is shown on the Figure [2] The chambers working in a single mode 
have a lower efficiency, which is well reproduced in the MC results. 
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Figure 2: A comparison between the input parametrization (solid line) and 
the MC results (dots). Left plot: efficiency distribution for all the RPCs; 
Right plot: example of the efficiency for all chambers from the fifth layer of 
one of the Barrel Wheels. Two of the chambers are working in a single gap 
mode, which is well reproduced by MC. 

Noise signals in the MC are defined as all the ffied strips without a 
corresponding simulated hit. For the purposes of validation, the average 
simulated noise for the chamber N^h has been calculated as: 
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where dnoise is the number of all noise signals for a given chamber in the 
event, Sch is the chamber the number of simulated events and 

T = Uev-ribx-'^^n-s is the total simulated time. 

The average simulated noise for the strip has been calculated in a similar 
way, but using the strip area and the number of the noise signals for a given 
strip. The comparison between the input noise parameters estimated from 
experimental data and the simulated noise is shown on the Figure [3j The 
time resolution, used in the simulation is 1.4 ns. Thus the timing is always 
correct, because it is measured in bunch crossing, where one bunch crossing 
is 25 ns. 




[Hz/cm2] RPC chamber 



Figure 3: A comparison between the input parametrization (solid line) and 
the MC results (dots). Left plot: Intrinsic noise distribution for all the RPCs 
strips; Right plot: comparison between the simulated noise and the input 
noise parameter for all chambers from the third layer of one of the Barrel 
Wheels. 

The comparison between cluster size distributions used as simulation 
input parameters and the MC results is presented on Fig. [4] (right) together 
with overall cluster size distribution (left). The cluster size depends on track 
impact position, relative to the strip edge. In our parametrization the strip 
is divided into 5 slices, i.e. the slice corresponds to the muon impact point. 

3.2 RPC system simulation vs experimental data 

RPCs have been mainly used as dedicated trigger detectors, however, they 
contribute to the muon reconstruction as well [21] . The main goal of the sim- 
ulation is to reproduce correctly the response of the RPC system, including 
efficiency, acceptance, misalignment, geometry, etc. To validate the correct- 
ness of the overall RPC simulation, the MC simulation has been compared 
to the real experimental data. The data is taken with the CMS detector 
during proton-proton collision runs with the energy y/s = 7 TeV at the 
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Figure 4: Cluster size for muons crossing RPCs. The mark presents the 
MC result and the line presents the input values. Overall cluster size is 
presented on the left plot. A comparison between input parameters (cluster 
size distribution) and MC result for five different muon impact positions 
(slices) is presented on the right plot. 

LHC in 2011. The environmental temperature during data taking is stable. 
The applied high voltage is corrected for environmental pressure and tem- 
perature variations |22] . The data analysis procedure is described in details 
in |21j . Here, the simulated average number of RFC hits associated to the 
muon has been compared to the number obtained from experimental data. 
The number of RFC hits used in the muon reconstruction is related to the 
number of RFC stations crossed by the muon (that is a function of rj) and of 
the intrinsic RFC efficiency. The comparison between the average number 
of reconstructed RFC hits along the muon track for data and MC is shown 
in Figure[5]as a function of rj and azimuthal angle ip. The shapes of the plots 
reflect the geometrical acceptance of the RFC muon system. The simulation 
model reproduces this feature well. 

4 Conclusion 

An algorithm for the simulation of the CMS RFC system has been developed 
and implemented in the CMS software package CMSSW. The set of simu- 
lation parameters has been defined using experimental data taken during 
cosmic and proton-proton collision runs. The parameter's values habe been 
regularly updated. A dedicated validation procedure has been developed. 
The MC and experimental data are in good agreement. 
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Figure 5: Left plot: The average number of RFC hits associated to the muon 
as function of ij for muons with Pt > 20 GeV/c from Z decay. Right plot: 
The average number of RFC hits associated to the muon as function of (p 
for muons with Px > 20 GeV/c from Z decay. Muons reconstructed both in 
the tracker and muon system are taken into account in this plot. The dots 
represent the experimental data and the solid line - the MC result. 
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